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Four-dimensional quantum field theories generally require regularization to be well defined. This 
can be done in various ways, but here we focus on Pauli-Villars (PV) regularization and apply it 
to nonperturbative calculations of bound states. The philosophy is to introduce enough PV fields 
to the Lagrangian to regulate the theory perturbatively, including preservation of symmetries, 
and assume that this is sufficient for the nonperturbative case. The numerical methods usually 
necessary for nonperturbative bound-state problems are then applied to a finite theory that has 
the original symmetries. The bound-state problem is formulated as a mass eigenvalue problem in 
terms of the light-front Hamiltonian. Applications to quantum electrodynamics are discussed. 
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1. Introduction 

In order to solve a (3 + l)-dimensional theory, the theory must be regulated in some way. In 
doing so, one should attempt to preserve as many symmetries as possible. We do this by adding 
enough Pauli-Villars (PV) [1] fields to regulate perturbation theory and assume that the nonper- 
turbative eigenproblem is also regulated. Numerical methods are then applied to a finite theory, 
just as was the case for (l+l)-dimensional superrenormalizable theories [2]. From the Hamiltonian 
eigenproblem, we can compute wave functions as coefficients in Fock-state expansions, and from 
these, compute observables. 

As a test of the approach, we consider QED. This is not meant to compete with perturbation 
theory; the numerical errors of the nonpertubative calculation are too large to resolve high-order 
contributions that perturbation theory computes directly. However, the method is intended for 
strong-coupling theories where perturbation theory is ineffective. 

There have been a number of applications of the method. The first being an exploration in 
terms of a soluble, heavy-source model [3, 4]. The dressed-fermion state in Yukawa theory has 
been studied extensively [5, 6, 7], as have exact solutions in the limit of equal PV masses [8]. 
Applications to gauge theories have been primarily to the dressed-electron state in QED [9, 10, 11, 
12, 13], but also to the photon eigenstate [14]. A scheme has been proposed for QCD [15]. 

In order to have a well-defined Fock-state expansion, the theories are quantized on the light 
front [16, 2]. We define light-cone coordinates of time, x + = t +z, and space, x = (x~ ,x±), with 
x" = t — z and x± = (x,y). The light-cone energy is p~ = E — p. and momentum, p = (p + ,p±), 
with p + = E + p- and p± = (p x ,p y ). These lead to the mass-shell condition, p 2 = m 2 , in the form 



We work with the standard parameterization, where the bare parameters of the Lagrangian are 
fixed by fits to physical constraints. There is the alternative of sector-dependent parameterization, 
where the bare parameters of the Lagrangian are allowed to depend on the Fock sector(s) on which 
the operators act. This was originally proposed by Perry, Harindranath, and Wilson [17] and applied 
to QED by Hiller and Brodsky [18]. More recent work with this approach has been by Karmanov, 
Mathiot, and Smirnov [19, 20]. Other nonperturbative approaches include lattice theory [21, 22], 
Dyson-Schwinger equations [23], effective -particle representations [24], and basis-function meth- 
ods [25]. 

The remainder of this paper contains a brief description of light-front QED in Lorentz gauge, 
followed by discussion of the eigenproblem for the electron state in lowest truncation. The anoma- 
lous magnetic moment is calculated. We also compare the results of the sector-dependent approach, 
and then conclude with a brief summary. 



2. Light-front QED in Lorentz gauge 

The Lorentz-gauge QED Lagrangian, regulated by two PV fermion flavors and two PV photon 
flavors, is 



= £(-iy 
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where 
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(2.2) 



A subscript of / = indicates a physical field, and i = 1 or 2 a PV field. The / = 1 fields are chosen 
to have negative norm. The mass of the bare photon /Xo is zero. 

The constants j3; and ^ control the coupling strengths of the various fields. We require that 
j3o = 1 and £o = 1 an d require the constraints L?=o( — 1)'A' = an d Tl=o(~ = 0- These guaran- 
tee the regularization and that the combinations i/a and in (2.2) have zero norm. A third pair of 
constraints comes from requiring that the photon eigenstate have zero mass [14] and that the mass 
of the electron eigenstate becomes zero when mo is set to zero [10]. 

The nondynamical fermion fields satisfy the following constraints (i = 0,1,2): 

j(-l) ; 'a_V^- + eoA-VftV- = iff 1 • [(-iyd±Wi+ -ie oy /piA ± Y+] - (-1)^7°^+- (2-3) 
From these we obtain a constraint on the null combination 

id-Y- = (iff) ■ d±Y+ ~ r°L^VftW'+- (2-4) 



The terms containing the photon field cancel because £,■(— 1 )'/$,- = 0; therefore, light-cone gauge is 
not necessary. The nondynamical field \j/ can then be constructed from a sum of that satisfy 
the free-fermion constraint. 

The mode expansion for the full Fermi field is 



1 



V167T 3 



dk 



b is (k)e 



-ihx. 
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(2.5) 



The spinors are defined in [26], and the nonzero anticommutators are 

{b is (k),bi,(k^)} = (-iys u ,s ss ,s(k-g) and {d is (k),dj sl (k^)} = i-iySffS^atk-k!). 

The mode expansion for the ith photon flavor is 



Vl6n 3 



dk r 



cn^e-^ + aJ^e 



ik-x 



(2.6) 



with the commutator [a^Qc) ,at (k^)] = {—\ydu!e ti d llv d{k — lc). The metric signature is chosen 
to be e 11 = (—1,1,1,1), for Gupta-Bleuler quantization [27]. Because we do not use light-cone 
gauge, there is no constraint on A + = A~ , and, consequently, there will be no instantaneous photon 
interaction term [2] in the Hamiltonian. 

We can now construct the light- front Hamiltonian from spinor matrix elements [14]: 



^ = E / dp rn ^(-V%(p)hs(p)+I t I dp^^{-i)%{p)dUp) 



2+. 



+Z[dkrf^(-l) l e%(k) ai ,(k) 
i,n J K 



(2.7) 
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+ I \JPMl I ^pdq{bl(p)[b J>s (q)V^(p,^ 

+ b J,-s(i) U U-2s(P> I) a \» (l ~ P) 

+ b l(p) [4M)Vil2s(pS+4-MWjMpS] M?+£) 

- d tsii) d jAi)vij&(p>i) + d jMq)u!}-2s(p,q) 4ii(i-p)+ Hx } ■ 

The vertex functions can be found in [14]. The Hamiltonian does not contain any instantaneous 
fermion terms [2]. They cancel between physical and PV contributions because they are indepen- 
dent of the fermion mass and proportional to ( — l)'j3; for the /th flavor. The sum over flavors then 
yields £,•(— l)'j3,' = 0- This is independent of the gauge choice. 

3. Electron eigenstate 

In the one-electron/one-photon truncation, the Fock-state expansion of the electron eigenstate, 
for total J z = ±i, is 



wHD) = E*4caio> + £ / dkc£ (k)bt(p-k)ai(k) |o) 



<JSfl 



(3.1) 



It is normalized according to {y a '(P!)\y a (P)) = 8(P — P)S a io- The second PV fermion flavor 
(i = 2) plays no role in this sector and can be removed. In order to have a positive norm, all 
physical quantities are computed from a projected state [9]: 



iVpUa> = E(-i)'«*L(p)i°>+E Ml I S 



(3.2) 



s\i - J (=0 y=0,2 

x I ( -^C^(k)bl(P-k)a^(k)\0)- 



To be an eigenstate of the light-cone Hamiltonian, the wave functions must satisfy the follow- 
ing coupled equations, with y = k + /P + : 



[M 2 -mj} Zi = j (P + ) 2 dyd 2 k ± £ ^(-iy + 'e^ 

x Vg CP - k, E)Cf ± (k) + Ufa (P - k, P)cg (k) 



(3.3) 



and 
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M 



m)+k\ tf + k 2 ± 



l-y 



eg® = v^I(-i)'"^p + ^ ± (p-_,p 



(3.4) 



(3.5) 



An index of i coiTesponds to the one-electron sector and j to the one-electron/one-photon sector. 
Therefore, in the sector-dependent approach, a mass m,- in a vertex function is assigned the bare 
mass, and ray is the physical mass. In the standard parameterization, all are bare masses. 
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The coupled equations can be solved analytically [9]. The wave functions C§ are 

nj+k 



«® = ^ M2 _w:w ■ (3 - 6) 



l-y y 

and the amplitudes satisfy 

(M 2 - m 2 )z ; = 2 e 2 £(-l) ! 'z i , [J+mimflo - 2(m i + m,,)/i] , (3.8) 



with [9] 



y " (Mj -y i6^ ^ m2 _^m_/^j(i-j)"' ( } 

l-y y 

j(m 2 ) = f^y Mil (3 10 ) 

l-y y 

There is also the identity [10] J = M 2 I . 
The analytic solution is [9] 

(M±m )(M±m l ) M±m Q 
ao± ~ 8^(mi -m )(2/i ±M/ ) ' 21 ~ Afimi 20 ' ( ' 

A graphical solution is given in [10]. In general, the lower sign yields the physical answer, because 
mo then becomes the physical mass M = m e at zero coupling. The coupling coefficient ^ is fixed 
by requiring that M = when mo = 0. In this truncation, we can safely take the m\ — > oo limit, 
where z\ = 0, mizi — > ±(M — mo)zo, and 

M(M±m ) 

ao± - ± 8^(2/ 1 ±M/o)' (3 " 12) 

and the second PV photon flavor can be discarded. In the sector-dependent approach, I\ and Tq are 
independent of mo, and the solution for Ofo can be rearranged as an explicit expression for mo 

m = TM + 8n < j^{2I 1 ±MI ). (3.13) 

The anomalous magnetic moment is computed from the spin-flip matrix element of the elec- 
tromagnetic current J + [28]. In the one -photon truncation and in the limit where the PV electron 
mass m\ is infinite, the expression for the anomalous moment is 

a e = —mlzl I y 2 {\ -y)dydk\ (3.14) 

K J 

X Q (-■)' V 

For the sector-dependent parameterization, the product (XqZq is just a, and the bare mass mo in the 
denominator is replaced by the physical mass m e . Plots of the anomalous moment as a function of 
the regulator mass jUi can be found in [12]. 
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4. Sector-dependent and standard parameterizations 

Even though we do not include the vacuum-polarization contribution to the dressed-electron 
state, the sector-dependent bare coupling is not equal to the physical coupling. Instead, they are re- 
lated by eo = e/zo [18], where zo is the amplitude for the bare-electron Fock state computed without 
projection onto the physical subspace. In general, the bare coupling would be eo = Z\e / y/Z^Z^/Z^; 
this includes the truncation effect that splits the usual Z2 into a product of different \f7^_ from each 
fermion leg [29]. With no fermion-antifermion loop, we have Z3 = 1, and without a second photon, 
there is no vertex correction and Z\ = 1 . Also, only the fermion leg with no photon spectator will 
be corrected by \[Z~2\ therefore, we find -sJZ^Z^f = Zo- 

In the sector-dependent approach, the bare-electron amplitude without projection is deter- 
mined, in the infinite-mi limit, by the normalization 1 = z\ -\-e\z\Ji, with 

j 2 = ydydk\Y(-\) k - ) , ' e . ^r . (4.1) 

y * X £ K 1 k 2 +l- y M +« 2 



Replacing eo by e/z.o, we can solve for zo as z.o = y 1 — e 2 J2 and find eo = ej \/\— e 2 Ji. For large 
jUi , one finds ^2 — ^1 + f ) • Thus, eo can become imaginary and Fock-sector probabilities 

range outside [0, 1] due to IR and UV divergences, and consistency then imposes limits on Ho and 
Hi [12], as confirmed for Yukawa theory in [20]. 

In the standard parameterization, the bare amplitude is determined by 1 = Zq + e 2 z\Ji, with 

h = 8^2 / ydydk\[ml-Amom e {\ -y) + m 2 (l -yf + k\] (4.2) 



x 



1 * 1 



ho ^ [kl + {l-y)H 2 k +ym 2 -y{\-y)m 2 ] 



Thus the bare amplitude is zo = 1/yl +e 2 J2, which is driven to zero as }X\ — > °° and causes most 
expectation values also to go to zero. Therefore, in this case there is a limit on pL\, but /lo can be 
zero. 

The anomalous moment in the sector-dependent case is 

a e = -m 2 [ y 2 (\-y)dydk\ (4.3) 

71 J 

1 



S, ( l)k {ym 2 + {\-y)ri+k 2 -m 2 y(l-y)) 



to \ym 2 + {l-y)lll+k\-m 2 y{\-y). 

In the jUi — > °°, no — y limit, this becomes exactly the Schwinger result 



a m2 r dyd^lQ-y) a ^ 



K 



III; 



In 



1-v 1 y e 

However, this limit cannot be taken without making the underlying theory inconsistent. 
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5. Summary 

With use of PV regularization, one can formulate and solve nonperturbative problems in field 
theories. It is important to maintain symmetries, which can be done with additional PV fields, 
such as those introduced to maintain a zero photon mass and the chiral symmetry of the massless- 
electron limit. It is best to regulate before applying numerical methods, to clearly separate limits 
of regulators from those of numerical convergence. The PV fields do add to the numerical load but 
also reduce it, by eliminating instantaneous fermion and instantaneous photon interactions. 

For both the standard parameterization and the sector-dependent parameterization, truncation 
of the Fock space results in uncancelled divergences, which require that not all PV masses be 
taken to infinity; however, meaningful results can be extracted at finite PV masses. For the sector- 
dependent approach, this is complicated by infrared divergences [12]. 

As discussed elsewhere [13], these methods have been extended to a truncation that includes 
two photons. The next step is to also include an electron-positron-pair contribution to the dressed- 
electron state and study charge renormalization as well as current covariance. This truncation will 
then include all contributions of order a 2 to the anomalous moment. A calculation at large a, where 
numerical errors in the order-a contribution would be small compared to the a 2 contributions, 
could be compared with higher-order perturbation theory. It would also be interesting to consider 
the dressed electron in a magnetic field and extract its induced magnetic moment. In addition, 
as a precursor to consideration of mesons in QCD, one can compute two-fermion bound states in 
Yukawa theory and QED. 
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